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Abstract. With recent advances in the precision of inclusive lepton-nuclear 
scattering experiments, it has become apparent that comparable improvements are 
needed in the accuracy of the theoretical analysis tools. In particular, when extracting 
parton distribution functions in the large-x region, it is crucial to correct the data 
for effects associated with the nonzero mass of the target. We present here a 
comprehensive review of these target mass corrections (TMG) to structure functions 
data, summarizing the relevant formulas for TMGs in electromagnetic and weak 
processes. We include a full analysis of both hadronic and partonic masses, and trace 
how these effects appear in the operator product expansion and the factorized parton 
model formalism, as well as their limitations when applied to data in the a; 1 limit. 
We evaluate the numerical effects of TMGs on various structure functions, and compare 
fits to data with and without these corrections. 
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1. Introduction 

The scattering of electrons on liadronic targets lias historically played an essential role in 
our understanding of the proton as a composite particle made up of partons: quarks and 
gluons [1-4] . Presently, data from electron and neutrino scattering at large momentum 
transfers, that is, deeply inelastic scattering (DIS), are used to determine the parton 
distribution functions (PDFs) which characterize the substructure of hadrons [5-10]. At 
lower energies, the resonant components of hadronic structure, and the duality between 
hadronic and partonic descriptions of interactions, continue to be explored [11, 12]. 

As the precision of the recent lepton-hadron scattering data has improved, it is 
vital for the theoretical analysis to keep pace. For example, the calculation of the 
Wilson coefficients has progressed to encompass next-to-leading order (NLO) quantum 
chromodynamics (QCD) and beyond [13-22]. It is important, therefore, to consider 
all sources of corrections which may contribute at a comparable magnitude, such as 
electroweak radiative corrections [23,24], quark mass effects [25-32], and target mass 
corrections [33-40]. In this review, we will focus on the problem of target mass 
corrections (TMCs), which formally are subleading 1/Q^ corrections to leading twist 
structure functions, where is the squared four-momentum transfer to the hadron. 

Understanding TMCs is important for several reasons. Their effects are most 
pronounced at large x and moderate Q^, which coincides with the region where parton 
distribution functions (PDFs) are not very well determined. A reliable extraction of 
PDFs here therefore demands an accurate description of the TMCs. Furthermore, a 
reliable interpretation of data on multiparton correlations at low momentum transfer 
depends on the proper accounting of TMCs. While target mass corrections have a 
long history, implementing these has not been entirely straightforward, as there exist a 
number of conventions, prescriptions, and potential scheme choices which can lead to 
differences in the final numerical results. 

The target mass corrections to electroweak structure functions were first determined 
by Georgi and Politzer in 1976 [34] within the operator product expansion (OPE) at 
the leading order of QCD. In the same year Barbieri et al. [35,36] rederived the mass 
corrections to scaling in DIS, including effects arising from non-zero quark masses. These 
same corrections were later derived from a parton model approach by Ellis, Furmanski 
and Petronzio [37,38]. Beyond leading order, the NLO QCD corrections to the target 
mass corrected structure functions were derived by De Riijula, Georgi and Politzer [39]. 
Recently, Kretzer and Reno [41, 42] reevaluated the TMCs for charged current (CC) and 
neutral current (NC) fi- and r- neutrino DIS, including NLO QCD corrections. 

There is a number of theoretical ingredients necessary to derive the TMCs to 
hadronic structure functions in lepton-hadron scattering in the context of the OPE, 
in order to relate them to the quark-parton model (QPM). The OPE method makes 
use of basic fundamental symmetries to relate the cross section to a reduced matrix 
element; as such, the OPE takes the hadron mass fully into account. In order to relate 
the reduced matrix elements to quantities which can be computed in the QPM, one then 
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inherits the associated hmitations. For example, the QPM describes the interaction as 
involving the scattering from a single, free parton (the leading twist contribution, see 
Sec. 2). Multi-parton correlations, formally higher twist, are discarded. Additionally, 
the QPM also imposes collinear kinematics on the parton involved in the interaction. 
The potential for the parton momentum to have a non-zero component transverse to 
the hadron momentum vector is neglected. This effect can omit mass contributions and 
introduce ambiguities if not addressed carefully. 

When using the QPM approach, the target mass corrections are taken into account 
in the following places: 

• in relating the parton fraction of the hadron's light-cone momentum (called the 
Nachtmann variable) to the Bjorken scaling variable x = Q"^ /2Mv, where M is the 
hadron mass and v the energy transfer; 

• in the mixing between the partonic and hadronic structure functions; 

• in a collinear expansion, where the TMCs appear from the px effects . 

In the limit of negligible target mass relative to Q^, the Nachtmann variable reduces 
to Bjorken-x. The QPM approach including TMCs has further limitations in the limit 
as X ^ 1. A "threshold problem" arises when trying to enforce that the structure 
functions vanish in kinematically forbidden regions. 

In the next section, we review structure functions and the OPE approach. We show 
the OPE results for the structure functions. Section 3 discusses the relation of the OPE 
to the parton model. We illustrate the formalism with the example of neutrino charged 
current scattering in Section 4. Section 5 describes the x ^ 1 problem and recent 
attempts to resolve this theoretical issue. Finally, in Section 6, we make some numerical 
comparisons using structure functions with and without target mass corrections. Our 
conclusions are summarized in Section 7. Three appendices detail the notation used 
in this review, the inclusion of charm mass corrections in neutrino charged current 
scattering, and a comparison of notation and results with Refs. [34,35,40] and [42]. 

2. Structure Functions and the Operator Product Expansion 

In this section we present an overview of the formalism for deep inelastic scattering 
(DIS). Using the framework of the operator product expansion (OPE), we outline the 
derivation of the leading electromagnetic and weak structure functions at finite Q^, 
including corrections from the non-zero target nucleon mass. We summarize these in a 
set of so-called "master equations" . Finally, we present an alternative but equivalent 
formulation of target mass correction (TMC) effects through the Nachtmann moments, 
and contrast this with the formulation in terms of Cornwall-Norton moments. 

2.1. Overview of Structure Functions and the OPE 

The basic lepton-nucleon inelastic scattering process, l{k) + N[P) i'{k') + X{Px), is 
shown schematically in Fig. 1, where k{k') is the initial (final) lepton four-momentum. 
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P is the target nucleon momentum, and Px is the momentum of the final hadronic 
state X. We define q = k — k' to he the four-momentum transferred from the lepton 
to the nucleon, with = —q^. The energies of the initial and final leptons are 
denoted by E and E\ respectively. Our notation reserves M for the target nucleon 
mass, P^ = M^, and the invariant mass squared of the final hadronic state is given by 
P^ = W^ = {P + qf = M^ + 2P-q-Q^ (see also Appendix A). 

For electromagnetic or weak neutral current (NC) scattering, the vector boson (V)- 
nucleon subprocess is V{q) + N{P) X{Px), where V = ■y, The related charged 
current (CC) process, which is important in neutrino-hadron scattering, i^(fc) + N{p) — »• 
i{k') + X{Px), where V = W^, will be discussed in Sec. 3, where we discuss the 
correspondence with the parton model. 

In addition to the virtuality of the exchanged boson, Q^, inelastic scattering is also 
characterized by the Bjorken scaling variable x, where 



In the massless target and quark limits (or equivalently in the Q"^ ^ oo limit), x is 
equivalent to the light-cone momentum fraction of the target carried by the interacting 
parton. In the target rest frame, the Bjorken variable can be written x = Q^/2Mi', 
where u = E — E' is the energy transferred to the hadronic system, and we define the 
inelasticity of the process hj y = u/E. For convenience, we also introduce the variable 
r to denote a frequently appearing combination of factors: 



At finite Q^, the effects of the target and quark masses modify the identification 
of the Bjorken x variable with the light-cone momentum fraction. For massless quarks, 
the parton light-cone fraction is given by the Nachtmann variable C, [33], 




(2) 



1 + y'l + 4x2MVg2 



(3) 
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X 

Figure 2. The Nachtmann variable ^ as a function of tlic Bjorken scaling variable x, 
for = 1, 2, 4 and 10 GeV^. For reference, a dotted line is shown for the limiting 
case = X. 



At large values oi Q^, C, ^ x. As Fig. 2 shows, however, for less than a few times the 
target mass of ~ 1 GeV, ^ can deviate significantly from x, especially at large x values. 
The Nachtmann variable appears naturally in the OPE, as we outline below. The full 
details of the notation, including parton masses, appear in Appendix A. 

We can write any generic inclusive lepton-nuclear scattering cross section as a 
combination of a hadronic tensor Wfj_i, and a leptonic tensor L'^'^: 

da ~ W^, L^' , 



where the hadronic tensor is given in terms of a product of hadronic currents,^ 

1 

^'-^w.^M^w.. (4) 

The structure functions Wi depend on x and Q^, as well as the target mass M. The 
hadronic tensor can be related to the discontinuity of the virtual forward Compton 
scattering amplitude T^j, via 

= -discT^, . (5) 

TT 

% In this work, we will focus on the unpolarized results. For TMC effects on the polarized structure 
function see Refs. [43-45], and references therein. 
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The OPE allows one to expand the hadronic matrix element in the forward 
scattering amplitude in a complete set of local operators [46]: 

T,,=zJ d'^z e*^-^ {N\T[J,{z) J.(0)]|iV) 

= T.<^;%-'"mN\OX...jN) , (6) 

i,T,n 

where the coefficient functions c^''^^,j"^"(q') represent the hard scattering of the boson 
from the parton. Here r denotes the twist of the operator O, defined to be the mass 
dimension minus the spin of the operator, and i represents different operators with the 
same twist. 

In the approximation of keeping only twist-2 operators, the forward scattering 
amplitude is explicitly: 

oo 
k=l 

22k 

X Qi^^2.n'^- '^- , (7) 

where A2k is the reduced matrix element of the twist-2 operator of spin 2k, and Cf' is 
the Wilson coefficient calculated using perturbative QCD. In Eq. (7), 

n»-. = v(-i).' '^"-{f ^ fci (py . (8) 

j=0 2^{2ky. j gt^nt^m's {2k~2j) pf^-^'s 

where {g.-.g} abbreviates a sum over (2/c)!/[2''j!(2A; — 2j)!] permutations of the 

indices. 

The strategy in the calculation, as outlined by Georgi and Politzer [34], is to 
evaluate Eq. (7), picking off each of the coefficients in the expansion, 

rp _ rp , PlJ.Pvrp ■ P'' Q rp , QfJ-Ql^ rp , P fJ-Qf ~^ PuQfl rp / n\ 

Finally, the structure functions Wt can be obtained from the imaginary parts of the 
amplitudes Tj using Eq. (5). In modern notation, the structure functions are denoted 
by Fi rather than Wi. To emphasize the inclusion of target mass corrections in the 
structure functions, we label Fi — s> F™'^, and note the correspondence 

f TMC 

{ Fii F2, F3, F4, F5 \ 

f 1 

\ 2xM^ xM^ 2M^ 2xM^ ''l " ^ ' 

In the following discussion, we focus on F™'~^ , F™^ and F™'~^ . The remaining two 
structure functions enter into the differential cross section suppressed by the lepton mass 
squared divided by ME [42], and can for most purposes be neglected. 

For a NLO calculation including an analysis of higher twist contributions sec. for example. Refs. [47- 

491 
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Using Eqs. (5) and (7), one can explicitly relate the Cornwall-Norton moments 
Ml'iQ'^) of the F™'^ structure functions to sums of reduced matrix elements [34,50]. 
For the F™'^ structure function, for example, one has: 



Jo 



k^K Q' J jKn - 2)! (n + 2j){n + 2j - 1) ' ^ ' 

Defining F^^^ to be the massless nucleon limits of the structure functions F-™^, we can 
relate the reduced matrix elements in Eq. (11) to the Cornwall-Norton moments of F}^'': 

= dyy-+'^-'F^'\y) , (12) 
Jo 

Cr^'A^^,, ^ f dyy-^'^-'Fi'\y) , z = 1, 3 . (13) 



In other words, the functions F^mc cQ^tain target mass effects, whereas the F^^^ do 
not. 

The Cornwall-Norton moments can be inverted to yield a set of target mass 
corrected structure functions derived from the operator product expansion: 

Q=) = f/iyo - '-^r^-^md)] , (14) 

Fr'=(^.Q^)=x^^{<i±^<„K)}. (15) 

F3™<^(x,0'^)=-.|-(i±^ft3(0). (16) 

The Fs^^C-expression comes from Eq. (4.17) of Ref. [34], and F™^ and F™^ follow 
in the same manner. The functions g2 and along with two more auxiliary functions 
hi^2i which will be needed below, are given by [42]: 

/^i(e,Q^)= /'rfn^^l!^^^, (17) 
h2{^,Q')= J^'du^-^^^^f^ , (18) 

hi^,Q')= Cdu ^' '^'^ ' , (19) 
92{^,Q') = j du h2{u,Q^) = j du I dv ^ ^ (20) 





u 






(m. 












(n. 





f2 Jf t)2 



Evaluating the derivatives in Eqs. (14)-(16) yields our final set of "master equations" 
for the target mass corrected DIS structure functions, which we discuss next. 
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2.2. Master Equations 

Combining the results in the previous section, the full, target mass corrected structure 
functions can be related to the massless limit functions by the following "master 
formula", using the notation of Kretzer & Reno (see Eq. (3.17) of Ref. [42], with their 
p r): 

i=l 

(21) 

where ^ is the Nachtmann scaling variable from Eq. (3) [33]. Inserting the coefficients 
Aj, Bj, Cj given in Tables I, II and III of Ref. [42], one finds for the first three structure 
functions: 



(23) 

F™'={x,Q') = I^Ff + ^-^h,(tQ') + , (24) 



with the functions hi{C,,Q'^) and g2i^,Q'^) given in Eqs. (17)-(20). The are the 



structure functions F™^ in the limit M — > 0: 

.(0) 



i^(°)(e,g^)= iimi^™^(x,g^ 



(25) 

=« 

Note that since ^ depends on x and M, F^^\^,Q^) ^ limM^o ^^'^(C, Q^), which has 
been the source of some confusion in the literature. Parton model representations of 
Fj'^^ will be shown in Sec. 4. 

We emphasize that the functions F™^ = F™'^{x,Q'^), and not F™^ = 
F™'^{^,Q'^), so that (x,(5^) is the correct point in phase space. While on the surface 
it may appear strange to have the left-hand-side of Eq. (21) be a function of x and 
the right-hand-side a function of ^, this arises quite naturally in the calculation. 
Specifically, evaluating the final state momentum conservation constraint, we can write 
(schematically) 6\q+P-Px) ~ S{x-0, and thus F™^{x,Q^) ~ Ff^\x,Q'^) 5{x-i) ~ 
Fi'^\^,Q'^)- Note that it would be incorrect to write F™^{^,Q^) ~ Ft\^,Q'^). All 
structure functions and PDFs depend on Q"^; we sometimes suppress this dependence 
for ease of notation. 

Another feature of Eq. (21) is that /12 and g2 appear in the formulas for both F™'^ 
and ^2'"'^*". This follows directly from the form of Eq. (7). For example, both the 
terms proportional to and to Cf'^ contribute to Ti (multiplying —g^y). The terms 
proportional to Cf'^' give rise to the second and third terms in Eq. (22). 

The "master equation" (21) holds to any order in the strong coupling constant a^, 
which implies that the coefficients A*, and Cj and the variable ^ are independent of 
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the order (LO, NLO, NNLO, . . . ) to which the structure functions -F/°^ are considered. 
In addition, Eq. (21) does not assume or imply any Callan-Gross relation. Specifically, 
one can compute the longitudinal structure function according to: 



4 



= + + ^vr^2(0 . (26) 

This general result gives a non-zero F™'^, and thus violates the Callan-Gross relation. 
The leading term (oc F^'^) is non-zero for finite quark masses, and the sub-leading terms 
(oc /i2,5'2) contribute for finite hadron mass M. 

Final state quark mass effects are taken into account by the parton model structure 
functions -F/'''*, and the general form of the master equation is unaltered. In other words, 
the {A*, Bj, Cj} coefficients and the Nachtmann variable in Eq. (21) will depend only 
on the hadronic mass and will not receive corrections due to final state quark masses. We 
illustrate this feature for the case of neutrino production of charm quarks in Sec. 4, and 
find agreement with results in the literature {of. Appendix C). In particular, the proper 
slow-rescaling variables (which depend on the quark masses) automatically appear as 
arguments of the parton distribution functions in a natural manner {of. Appendix A). 
For the general case of non-vanishing initial and final state quark masses we assume 
the same pattern holds true — the quark masses appear within the and the form 
of the master equation remains unchanged; this is a consequence of factorization. This 
ensures the correct massive parton model expressions are, by construction, recovered in 
the M — > limit; a necessary condition for any formalism which includes target mass 
effects. The modular structure of the master equation renders computations of target 
mass corrections simple and transparent once the parton model expressions for the 
(with or without quark masses, in leading or higher order) are given. 



2.3. Nachtmann Moments 

An alternative, but closely related, formulation of the TMCs is in terms of Nachtmann 
moments. Nachtmann [33] showed that one could arrange the OPE so as to ensure that 
at a given order in only operators of a given twist would appear. The Nachtmann 
moments /i^ of structure functions F™'^ [i = 1,2,3) are constructed from operators of 
definite spin. This means that from the infinite set of operators of twist-two and different 
spin contained in the trace terms of the OPE, only the operators of spin n contribute 
for the n — 2 Nachtmann moment of the structure function. This is contrasted with the 
Cornwall-Norton moments in which different spin operators contribute to the twist-two 
moment. The Nachtmann moments are defined to factor out the target mass dependence 
of the structure functions in a way such that its (Cornwall-Norton) moments would equal 
the moments of the corresponding parton distributions. 
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For the F^mc structure function for example, the Nachtmann moment is 
constructed so that it depends only on the reduced matrix element An (and Wilson 
coefficient C") on the right-hand-side of Eq. (11), in contrast to the infinite series of 
{M'^/Q'^y terms in Eq. (11) for the Cornwall-Norton moment. In effect, the subleading 
{M'^/Q^y terms with j > 1 are absorbed in the redefined moment on the left-hand- 
side of Eq. (11), so that only the first term with j = contributes to the Nachtmann 
moment. 

Specifically, the Nachtmann moment of the F2 structure function is given by: 



= / dx- 

Jo 



(n + 2)(n + 3) 

where again we use r = ^Jl + Ax'^M'^ /Q"^. One can also express the Nachtmann moment 
as an integral over the Nachtmann variable 

r«o ^„„2 (1 + eMVQ'f 



Jo 1 



3(r- 1) _ 3(r- 1)^ 
r2(n + 2) ~ r2(n + 3) 



^2M2/Q2 

X F™^{x,Q'^) |x=c/(i-Af2c2/Q2)' (28) 

with F™^{x, given by Eq. (23), and = ^{x = 1) = 2/(1 + y^l + 4M2/Q2). 

Similarly, for the longitudinal Nachtmann moments (or equivalently for the F™'~" 
structure function, defined in terms of F™'^ and F™'^), one has 

fi ^^^^ 



f^UQ')= j dx^-^\Fl^^{x,Q-) 



+ (n + 2)(n + 3) i^^) 

Note that while in the Q"^ ^ 00 limit fi^lQ"^) approaches the Cornwall-Norton moment 
of F^^c, at finite both F];^^ and ^2^^^^ contribute. 

The Nachtmann moment /i2(Q^) ^ind Cornwall-Norton moment M^((5^) can be 
related by expanding the moments in powers of 1/Q^. Expanding to 0{1/Q^), one 
has: 

n(n^-l) M^ _ n{n^-l) M^ 

+ 2(n + 3) Q4^^^2 W ) g + , (dUj 

which illustrates the mixing between the lower and higher Cornwall-Norton moments. 

In Fig. 3 we compare the n = 2 Nachtmann moment /i2^^ of the target mass 
corrected proton F™'^ structure function (solid curve) with the Cornwall-Norton 
moment M^=2 (dotted curve) of the same structure function in the massless target limit, 
^2°^ (note the expanded vertical scale!). While the two moments agree well at large Q^, 
a clear deviation from equality is seen for <^2 GeV^. Part of this discrepancy may be 
attributed to the behavior of the target mass corrected structure function in the x — 1 
limit, as we discuss in detail in Sec. 5. 
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0,205 




Q' (GeV") 

Figure 3. Comparison of the n = 2 Cornwall-Norton moment of the proton ^2"'' 
structure function (solid), with the Nachtmann moment of the target mass corrected 
F2™^ (dotted), calculated from F^°'' using Eq. (21). 



3. Relation of the OPE to the Parton Model 

In the previous section we used the OPE to relate the inclusive lepton-nuclear scattering 
cross section to the structure functions F^'^^ of Eq. (21). While these relations are quite 
general, they are of no utihty unless we can relate them to calculable quantities. Up to 
this point we have not invoked the parton model and its associated assumptions; the 
only ingredient has been the OPE, which only makes use of fundamental symmetries. 
In this section we briefly discuss the relation of the leading twist OPE treatment of 
inclusive lepton-nuclear scattering with the parton model. We also comment on the 
importance of TMCs in testing the validity of leading twist descriptions of data, or 
alternatively, extracting higher twist contributions. 

The leading, j = 0, term in Eq. (8) reduces the target mass corrected structure 
function F™'-^ to the massless limit function -F/^"*, which can be expressed in terms 
of parton distributions in the parton model. The j > terms in Eq. (8) resum the 
target mass corrections. The explicit form of the structure functions f/*^^ depends on 
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the interaction. For Fi'\ for example, the j = term in Eq. (13) can be schematically 
written (neglecting quark electroweak charges) as: 

= [' dxx''-'Fi''\x,Q^) (31) 
Jo 

= [' dzz--'a{z, Q\ /i^) . rfyy"- V(y, /i') (32) 

JO JO 

= dxx-' C ^a(-, Q\ ^^') /(y, ^^') , (33) 
JO Jx y ' 

where / is the parton distribution function, defined at a scale /i, and a is the boson- 

parton scattering cross section. The first factor in the product C" ■ An above is the 

process-dependent Wilson coefficient, while the second is the moment of the parton 

distribution /. At leading order in QCD perturbation theory, the parton cross section 

a{z,Q'^, fi"^) ~ (5(1 — z), multiplied by numerical factors that account for the process 

(such as e'j/2 for electromagnetic scattering with a parton of electric charge Cf). QCD 

corrections introduce a dependence to a and therefore to C". Both and A"' 

depend on fi^ in such a way that the product is /i^ independent. 

As remarked above, Eq. (21) is valid at leading twist, and to all orders in as- The 
inclusion of higher orders in f/^"* is standard in the leading twist approximation. This 
involves evaluating the perturbative QCD corrections to Vq q and associated gluon 
processes. The NLO corrections have been known for some time [13-16,27], and the 
NNLO corrections have also been evaluated [17-19,21,22,51]. 

In addition to the perturbative corrections, moments of structure functions in the 
OPE receive l/Q"^ power corrections, which are suppressed at large Q^, but may be 
significant at lower Q^. The size of these corrections is determined by matrix elements 
of local operators which have twist greater than 2, denoted "dynamical higher twist". 
Since one of the motivations for studying TMCs is to more reliably extract leading 
twist parton distributions from structure function data, it is necessary to ensure that 
TMCs are properly taken into account so as to reveal genuine effects associated with 
dynamical higher twists. Of course, the TMCs themselves enter as 1/Q^ corrections, 
despite formally being of leading twist, so in practice it is crucial to remove from the 
data TMC effects which could otherwise resemble higher twists. 

The general expansion for the moments of structure functions in QCD is then one 
which involves an expansion in powers of l/Q"^, each term of which has an associated 
perturbative expansion, in addition to the {M'^/Q'^y TMC expansion as in Eq. (11). 
For the above example of the Fi structure function, the general OPE expansion in the 
M ^ limit can be written [46]: 

I dxx-'Fi'\x,Q') = c-mAn + E^r(Q^) f=j 



n,i 1 



(34) 

where the first (twist-two) term corresponds to that in Eq. (31), and the higher twist 
coefficients B^i are of the same order as A^, with Mq some typical hadronic mass scale. 
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Figure 4. An example of a higher twist diagram, involving a four-quark operator, 
contributing to the term Bnj in Eq. (34). 

Several analyses of low data have been made which include both leading and higher 
twist contributions, although the practical difficulty in disentangling l/Q"^ from even 
higher order corrections means that usually terms only up to twist 4 are considered. 

Although in the context of parton distribution analyses higher twists are seen as 
unwelcome complications, the study of higher twists is also important in its own right. 
Because the higher twist operators necessarily involve several quark or quark and gluon 
fields, their matrix elements can provide insight into nonperturbative, multi-parton 
correlations in the nucleon, and possibly even on the long distance partonic interactions 
associated with confinement. For example, twist-4 corrections to the F2 structure 
function involve the so-called "cat's ears" diagrams schematically shown in Fig. 4, which 
represent flavor non-diagonal transitions between the incoming and outgoing photons in 
the Compton scattering process [52]. Differences in their magnitude in the proton and 
neutron can reveal the mechanism responsible for the phenomenon of Bloom-Gilman 
duality in structure functions [12]. For spin dependent structure functions, certain twist- 
3 and twist-4 matrix elements can be related to color polarizabilities of the nucleon, 
which describe how the color electric and magnetic gluon fields respond to the spin of 
the nucleon [53-56] . Before we can begin to unravel these effects, it is essential that the 
kinematical TMCs are properly removed from the data. 

Having outlined the formalism, in the next section we display explicit expressions 
for for leading order QCD at leading twist for neutrino-proton scattering. We focus 
on this process because of the subtleties associated with charm production by massless 
quarks. Next-to-leading order expressions for the structure functions in a fixed 
flavor number scheme with three active flavors (3-FFNS) can be found in Ref. [41]. 

4. Corrections for Finite Hadron and Quark Masses 

To illustrate how the OPE-parton model formalism is implemented in the case of a 
quark mass scale ruj in the flnal state (in addition to the target mass scale M), we 
consider the case of neutrino inclusive lepton-nuclear scattering for 3 light {u, d, s} 
flavors and one heavy {c} flavor. The reaction takes place via the process u^N fi~X 
(or iJfj_N /i+X), of. Fig. 5. In this example we will focus on the leading order only. 
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1^" _ 1^ 




Figure 5. The basic Fcynman diagrams for the DIS process for neutrinos, ^ X 

(left diagram) and anti- neutrinos DN —f fj.'^X (right diagram). 

SO we can separately observe how both M and rrij = rric enter; once the pattern is 
developed here, the generalization to multiple quark masses {mi,mj) will be outlined in 
Appendix A. 



4-.1. Multiple Mass Scales in Neutrino DIS 



To illustrate how the quark and hadron masses enter into the structure functions, we 
separate the M — > limit structure functions into contributions with a light (/) and a 
heavy quark (c) in the final state. For neutrino scattering, we thus have: 

F-'^'\^,Q') = F':f\i,Q')+Flf{i,Q') . (35) 

For neutrino DIS with a light quark in the final state (ie., no charm), the structure 
functions for a 3-flavor number scheme are given by: 



3,/ 



+ smvus?+u{myud? + \Vus?) 



dmv^d? + s(oiKsi' - m{\Vud? + IK, 



(36) 
(37) 
(38) 



where Vij are the CKM matrix elements, and the dependence in the PDFs has been 
suppressed. For heavy quarks, the structure function contributions are: 

.;.{0) 



(39) 
(40) 
(41) 



For anti-neutrino structure functions F^'^^\^,Q'^), the results are obtained by 
interchanging q ^ q in Eq. (35): 



-^1,2 



F^fiq^q], F. 



P,{0) 



F3'-°\q^q], q = u,d,s 



(42) 



The variable ^ in Eqs. (39)-(41) is the so-called slow rescaling variable, generalized 
to include both target and quark mass effects, and corresponds to the light-cone 
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fractional momentum of a massless quark which produces a charm quark. It is related 
to the Nachtmann scaling variable ,^ by = ^Rij, where Rij contains the dependence 
on the initial and final quark masses {of. Eq. (A. 8)), and in the limit of a single heavy 
quark j = {c}, and three light flavors i = {u,d,s} in the initial state, it reduces to 
Rij = (1 + ml/Q"^). In this case we recover the relation: ^ = ^ (1 + m'l/Q'^) [25, 29, 57]. 

Note that Eqs. (36) and (37) lead to a generalized Callan-Gross relation at leading- 
order (with X replaced by 2^F^'i'^\^) = F2f^\^). It is also important to note in 
Eqs. (39), (40), and (41) that operationally F^^'^^ is a function of the variable ^, while 
the PDFs are functions of = ^Rij. As discussed in Appendix B, this is simply a 
consequence of evaluating the final state partonic momentum conserving delta-function 
6{^-^/Ri,). 

4-2. Relation to TMC Structure Functions 

From the master equation (21), we can relate the above partonic distributions to the 
TMC structure functions. The charm mass effects can be taken into account in the 
structure functions F^^^'' in Eqs. (39)-(41) by using the generalized Nachtmann variable 
^, thereby allowing the form of Eq. (21) to be kept unmodified. For example, in the limit 
of a diagonal CKM matrix, Vij = 6ij, and neglecting the non- leading terms, inserting 
F2c°\^j Q^) — 2^s(,^) into the leading term of Eq. (23) gives: 



3 ■ </Qr . (43) 



The complete set of expressions in this "diagonal CKM" limit for the charm 
production component of the target mass corrected structure functions is given by: 

F™^(x, Q') ^ f (1 + mimsiO + ^(1 + f_ du ^ , 

4MV /•! /-i s{v) 

an 



24MV /•! /-i s{v) ,^ , 

' du dv , (45) 



Q^r^ Je, Ju V 

i^jf "(X, Q') - ^(1 + mimno + ^ Cdu^-^+o. (46) 

t,r^ Q^r-^ u 

A discussion of the dependence of hi and g2 on the charm quark mass is given in 
Appendix B. * 

Having illustrated the implementation for non-zero hadron and quark masses, the 
extension to the general mass case is straightforward from the relations in Appendix B. 

* See also Rcf. [42], Eqs. (3.1)-(3.3); with g2 and /13 defined as in Eqs. (B.4) and (B.5), the charm 
mass dependent structure functions of Eqs. (44)-(46) are recovered. 
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As with the above charm mass case, the form of the equations will remain the same; the 
target mass dependence enters in the replacement x — > ,^ = xRm, where Rm = 2/(1 + r) 
(see Appendix A), and the quark masses enter via the replacement ^ — > ^ = C,Rij- 

5. Threshold Effects and the x ^ 1 Limit 

The standard TMC formulation described in the previous sections involves the so- 
called "threshold problem", which was recognized soon after the original derivation 
of TMCs, and has been discussed by a number of authors in the literature [39, 58- 
61]. It is associated with the behavior of parton distributions in the threshold region, 
between pion-pro duct ion [W = M + m-j^) and the elastic point [W = M), and becomes 
increasingly important as x ^ 1. 

The problem can be summarized as follows. For simplicity, we neglect perturbative 
QCD (pQCD) corrections, with the Wilson coefficient functions Cf set to unity. 
Following Ref. [34], we define a leading twist parton distribution in this discussion by 
F(^,Q2) = Ff^(^,Q2)/^2 ^^2)). The n-th moment of the parton distribution, 

A^= f'd^CF{^,Q') , (47) 
Jo 

should be independent. Neglecting higher twist (HT) contributions, one should have 
[39,58]: 

/' (e, Ql) = t di C~' i^f (e, Ql) [no pQCD, no HT] (48) 

Jo Jo 

for any two momentum scales and Q2, where F2'\^, Q"^) is the structure function in 
the massless target limit, M — s> (see Eq. (25)). Since F!f'\^^Q'^) must vanish in the 
kinematically forbidden region ^ > ^o? where 

UQ^) ^ax = l, Q') = 2/(1 + v/l + 4MVg2) , (49) 

the equality in Eq. (48) implies that -p2°'*(^, Q^) must be zero for both ^ > ^o(Qi) and 
^ > ^oiQl)- If Ql < Ql in which case ^o(Q?) < ^o(<52)' this implies that Fi^\^,Ql) 
must vanish in the range ^o(Qi) < ^ < ^oiQl)- However, there is no physical reason in 
the parton model for it to do so here, and this leads to an unphysical constraint. 

According to the prescription of De Rujula, Georgi and Politzer (DGP) [39], the 
solution to this "paradox" lies in the higher order terms in the twist expansion for 
the moments of -F2°^(^,Q^). In general, the n-th moment of the structure function 
F2°\C,Q^) can be written as [39] {of. Eq. (34)) 

I d^ i^f (e, Q') = + E l^^j B^, , (50) 

where the i-th term in the sum is the contribution of operators of twist 2i + 2. 

DGP argue [39] that the higher twist B^^i terms indeed tame the paradox. They 
note that if the distribution at large ^ behaves as (1 — S,Y, with a ~ 3 — 4, then 
the n-th moment of f!^^ will be sensitive to the distribution at values of ^ close to 
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^„ = n/(n + a) ~ 1 — a/n for large n. For the moments to be sensitive to the threshold 
region, one must take n large enough so that C,n ~ Co(Q^) ~ 1 — M^/Q"^ at large Q^. 
In other words, the analysis is affected strongly by the region in moment space where 
n ^ aQ'^/M'^. For such large values of n, however, the terms proportional to Bn^i cannot 
be disregarded, and DGP conclude that there is therefore no paradox. 

More specifically, there is a non- uniformity in the limits n ^ oo (or ,^ ^ 1) and 
oo. The appearance of powers of nM^/Q"^ in Eq. (50) signals the breakdown of 
the whole twist expansion in this region. The master equation (21) is therefore not valid 
for large ^, where higher twists are important. 

The problem with the nonuniform limits was also recognized by Tung and 
collaborators [59,60], in the context of Nachtmann moments, which were discussed 
in Sec. 2. In particular, they pointed out that the Nachtmann moments n^{Q'^) of 
structure functions have, for fixed n, the asymptotic — > oo limit (again, neglecting 
pQCD corrections) 

/i^(g2) ^ An [n fixed, oo, no pQCD] . (51) 

However, the asymptotic form does not exhibit the correct threshold behavior at fixed 
Q^, which should be [59,60] 

/^r(Q') - QiQ') J^nQ') IQ' fixed, n^oc], (52) 

where Jx^{Q'^) are the "regularized" moments of a function which has support over the 
correct physical range, < < .^o- In the asymptotic n ^ oo limit, the regularized 
moments are then weakly dependent upon n, while the full moments /i^ contain the 
main n dependence through the threshold factor C.^iQ'^) [59,60]. 

Tung et al. [59, 60] proposed a solution to this problem by postulating an ansatz 
for the moments in which the threshold factor is explicitly included in the definition, 

/ir^^"^^^(g^) = C(Q') A., (53) 

where the new moments are consistent with both the asymptotic QCD behavior (51) and 
the kinematic threshold requirement (52). However, as they note, such a prescription 
is not unique, and in fact agrees with the standard OPE expansion only in the ra — oo 
limit [61]. 

The region of where TMCs are significant also corresponds to the region where 
nucleon resonances play an important role, W <^ 2 GeV. Although the resonance 
region displays significant dependent structure, a twist expansion is still useful here 
when one averages over individual resonances. The resonance region corresponds to 
n ^ aQ'^/M'^, and here higher twists are small but not negligible. The size of the higher 
twists in fact determines the degree to which the physical structure function oscillates 
around the leading twist function, in the sense of local Bloom-Gilman duality [11]. 
Target mass corrections are therefore closely related to the physics of resonances and 
quark-hadron duality [12]. 

Duality is also invoked implicitly in the DGP approach in that the parton 
distribution in the region between the pion production threshold, corresponding to 
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W = M + m^, and ^ = 1 is taken to be dual to the elastic form factor at ^ = 
[11, 62, 63]. However, a consequence of a non-zero parton distribution F{C,, Q^) over the 
entire region < ^ < 1, including the unphysical region ^ > (at x > 1), is that the 
target mass corrected structure function is non-zero at x = 1, i.e., F™^{x = 1, Q^) > 0, 
for any finite Q^. This seems necessary in the DGP approach if one wishes to preserve 
a probabilistic interpretation for F{^), as in the parton model. 

An alternative to working with distributions in the unphysical region > .^o w^^s 
suggested by Steffens and Melnitchouk (SM) in Ref. [61]. The philosophy adopted there 
was to reformulate the DGP analysis in such a way as to ensure the correct kinematic 
limit for the target mass corrected structure function, i.e., F™'^{x 1,Q^) 0, 
for any finite Q^. The method specifically involved working with parton distributions 
which vanish beyond the kinematic upper limit, ^ = C,o, and which at finite therefore 
depend on ^ and ^o- 

The approach of SM [61] defines the matrix element A^^^ in terms of an integral 
of the parton distribution over the range < < ^o? 



where the function F is now ^ and .^o dependent, and vanishes smoothly as ^ .^o? ^oi 
instance as a power of (■Co ~0- The function F is also properly normalized; for example, 
if the functional form for F is F{^, ^o) = N C,"" — then the normalization constant 



where now the function F depends only on a single variable ^. In this approach, 
one does not need to invoke higher twists to cancel the leading twist TMCs in the 
^ > ^0 region. In the original DGP formulation, the extended, unphysical range ^ > 
implies that the target mass corrected structure function is overestimated at large ^. A 
resulting complication of this is that even neglecting pQCD, the moments A^^^ become 
dependent, thereby making the connection with the OPE, and hence a partonic 
interpretation, problematic in the presence of TMCs. This dependence is, however, 
in accordance with the results of Tung et al, Eq. (53). 

An alternative approach to avoid introducing this additional Q'^ dependence is to 
redefine the moments such that 




(54) 



is AT = eo"""^"'. 



One can also define a "scaled Nachtmann variable' 




(55) 



(56) 




(57) 



In terms of the variable ^ this can then be rewritten as 




(58) 
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Figure 6. Ratio of the n = 2 Nachtmann moment of the F2 structure function and the 
n = 2 moment of the quark distribution, as a fimction of Q^. The curves correspond 
to standard TMC prescription [34] (dotted) and the modified TMC (SM) prescription 
from Ref. [61]. 

This definition leads to the standard approach of Georgi & Pohtzer, A^^''^ An, 
with the unphysical region contributing in the calculation of the physical structure 
functions. Hence, we arrive at a conundrum: imposing the correct threshold leads 
to a dependent An, in contradiction to the standard OPE; retaining the partonic 
interpretation of An gives rise, on the other hand, to unphysical contributions to the 
physical structure functions, which one hopes will be canceled by higher twists. 

The large-x behavior of the distributions is also particularly important for the 
comparison of the Nachtmann moments of structure functions with the moments of 
the parton distributions. The Nachtmann moments are, by construction, meant to be 
protected from target mass effects, so that these should be equal to the moments of 
the scaling parton distributions for any Q^. A comparison of the Nachtmann moment 
with An was already shown in Fig. 3, and in Ref. [61] this property was investigated 
in several other scenarios including the modified definition of the parton distribution in 
Eq. (56). 

In the standard approach, Eq. (58), there is an increasing difference between the 
Nachtmann moments and the moments of the parton distributions for < 1 GeV^, 
as can be seen in the ratio shown in Fig. 6. At Q"^ = 1 GeV^, the difference is ~ 2%. 
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In the modified approach, on the other hand (Eq. (56)), the Nachtmann moments are 
essentially equivalent to the .^o-dependent A^^^ for as low as 0.5 GeV^, as shown 
with the solid line. 

The deviation of the "standard TMC" ratio from unity below ~ 2 — 3 GeV^ 
suggests that the equality between the Cornwall-Norton moments of the PDFs (or 
massless limit structure functions) and the Nachtmann moments of the target mass 
corrected structure functions will only hold at sufficiently large Q^. For > 2 GeV^, 
we can expect uncertainties due to the treatment of TMCs to be 1%. These effects will 
also appear when the evaluations are done using experimental data. For lower values of 
Q^, the direct connection between Nachtmann moments of F™'^ and the parton model 
requires a more sophisticated treatment of the large-x regime. 

6. Quantitative Effects of Target Mass Corrections 

The target mass corrections described in the previous sections, and summarized in the 
master equations (21), are applied in this section to neutrino-nucleon and electron- 
proton scattering. The theoretical evaluation of the target mass corrected structure 
functions F2 and F3 are compared with data from the NuTeV neutrino-iron scattering 
experiment [64] . (In this discussion we neglect nuclear effects on structure functions, and 
assume equality between the nuclear and nucleon structure functions — for a discussion 
of nuclear effects in inclusive lepton-nuclear scattering, see e.g. Ref. [65].) We also show 
results for the electromagnetic structure function in electron-proton scattering, with 
a comparison data. We begin in Sec. 6.1. with the ratio of F2 and F^ evaluated in uN 
scattering {N = {p + n)/2) with and without target mass corrections. 

6.1. TMC Effects in the Massless Quark Limit 

Target mass corrections are relevant especially in the high-x and low-Q^ regions. To 
quantify the size of the corrections, we evaluate the structure functions at leading order 
in QCD, for neutrino-nucleon scattering in the massless quark limit. For convenience 
we use the CTEQ5 PDF set [6]; using other parameterizations has a minimal effect on 
our numerical results. The effects of the target mass corrections are illustrated in Fig. 7, 
where the ratios of the F2 and F^ structure functions with and without TMC effects are 
shown, namely F™^/fP, z = 2, 3, for = 1, 5, 25 and 125 GeVl 

The ratios of structure functions with TMCs to those without rise above unity at 
large x, with the rise beginning at larger values of x as increases. The correction can 
be quite large: for x = 0.8, for example, the TMC effect is ~ 30% at g2 = 5 GeVl A 
large part of the correction comes from the shift of x ^ < a; in the PDFs. The PDFs 
are rapidly falling functions at large x, so even a small change in the argument of the 
PDFs {x replaced by C,) can have a sig nificant impact on F™^ / f'>^\ For Q"^ = 1 GeV^ 
^ deviates from x even for x ~ 0.3 (see Fig. 2). In the following section we discuss the 
relative importance of the non- leading terms in the master equations (21). 





6.2. Non-leading Terms in the Master Equations 



The non-leading (second and third) terms in Eqs. (22)-(24) constitute a small correction 
to the leading term. Since the evaluation of the convolution integrals is quite time- 
consuming, it is useful to have an upper bound for the size of these terms. We can 
rewrite Eq. (23) as: 



X 



1 + 



with n = M'^/Q'^. The structure function ^2°^ appearing in the integrals in Eqs. (18) 
and (20) is a decreasing function of x or ^ (see e.g. Fig. 10 below). Consequently, 



F2°'* can be evaluated at the lower integral limit, giving /i2(0 < {^2^' iO/Oi^ " 
92(0 < -^2°^(0(~ I'^'C ~ 1 + 0- O'^^ then arrives at the following inequality: 



.(0), 



?TMCi 



X 



x,g^)<-^Ff (0 



1 + 



Ine 



1 + 

(59) 
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The expressions (6/ix^/r)(l— ^) and (12/i^x^^^/r^)(— In^ — 1+^) can be easily evaluated 
to obtain an upper bound for the contribution of the non-leading terms. Following the 
same line of argumentation one finds for the structure function F3: 



(60) 



While the upper bounds for ^2™^ and F™*" are strictly satisfied for x and values 
relevant for target mass corrections, these bounds are of limited practical use. For 
example, for Q'^ = 1 GeV^, Eq. (59) places a limit on the non- leading corrections to be 
less than ~ 65% of the leading term at large x. The actual value is much less, below 
~ 21%. Therefore it is useful to note that 



„2 



6fix^ 



(61) 



provides a very good approximation of the structure function F™'"{x,Q'^). Similarly 
F™'^{x,Q'^) can be approximated by 

fixC, 



F3™^(x,Q^)^|,Ff(0 



-(1-Oln^ 



(62) 



The magnitude of the non-leading contributions to the target mass correction is 
illustrated in Fig. 8, where the ratio of the target mass corrected F2 (top graph) and 
-F3 (bottom graph) structure functions is shown relative to the leading contribution, 
at Q2 = 1^4 and 10 GeV^. Here F2™'^'"^(x, Q^) = (xVe'^')^2°^(e) represents 
the leading contribution to the target mass corrected structure function F™'^{x,Q'^) 
{cf. Eq. (23)), while the corresponding results for the F3 structure function is given 
by F™'^'^'^^'^^^^{x,Q'^) = {x / C,r'^) F^'^\c,) . For definiteness, the structure functions are 
for charged current neutrino-proton scattering and have been computed in next-to- 
leading order of QCD including quark mass effects. However, the results are very robust 
concerning variations of these details (process, order, quark mass effects). The dotted 
curves in Fig. 8 present the results of the approximate formulas in Eqs. (61) and (62). 
As can be seen, the simple approximations are in very good agreement with the exact 
results. 

The excess over unity depicts the fractional contribution of the non-leading terms. 
Clearly, the non-leading contributions to the structure function F™'~' are relatively small 
and positive. For = 1 GeV^ they amount about 21%. However, for = 4 GeV^ 
they correct the leading term already by less than 6%, and for = 10 GeV^ by 
less than 3%. The major contribution from the non-leading pieces comes from the 
terms proportional to /i2 or h^, whereas the part proportional to (72 constitutes a small 
correction. These results imply that F™'^{x,Q'^) can be approximated in many cases 
by F™'^'^^^'^^^^^{x, Q^). Moreover, if more precision is needed, the non-leading pieces can 
be simulated, to good accuracy, by the simple approximation in Eq. (61). For F™'^, 
the non-leading terms are even smaller, contributing less than 7% for = 1 GeV^. For 
= 4 GeV^ they are already below the 2% level. 
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Figure 8. Ratio of the target mass corrected F2™^(x, Q^) (top) and i^3™°(x, Q^) 
(bottom) structure functions to the leading contributions in Eqs. (23) and (24), at 
= 1,4 and 10 GeV^. The soHd curves represent the exact results, while the dotted 
curves have been calculated using the approximate formulas in Eqs. (61) and (62). 



The findings of this section are useful in providing a simple explanation of the main 
effects of the target mass corrections shown in Fig. 7. In fact, to good accuracy, the 
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ratio in Fig. 7 is given by: 

Since ^ < x, and F2 is a monotonically decreasing function, one can see that the simple 
shift from x to ^ in the argument of ^2°'' (along with the factor x^/^^ > 1) produces 
the enhancement over unity visible in this figure. The dip below unity for low values 
and intermediate x can be explained by the factor in Eq. (63) which is smaller than 
one. 



6.3. TMC Effects in NuTeV Structure Functions 

A comparison of the F2 structure function with the NuTeV neutrino data [64] is shown 
in Fig. 9, where F2 is calculated both with (solid curves) and without (dashed curves) 
target mass corrections. For this comparison, the full NLO evaluation is made including 
charm quark mass effects. The PDF set for this comparison was taken to be CTEQ6HQ 
[7]. As Fig. 9 shows, there is a clear improvement in the agreement in the high-x region 
between the NuTeV data and NLO predictions including target mass effects. 



6.4. Unfolding Target Mass Effects From Structure Function Data 

In neutrino and charged lepton scattering from nucleons it is the F™'~^{x, Q^) functions 
that are measured by experiment. In principle, measurements over a range in x at 
fixed allow for a determination of the massless limit structure functions Fi'^\x) by 
inverting the master equations of Sec. 2. Implicit in such a procedure is the assumption 
that dynamical higher twists operators either contribute negligibly to the data in the 
region of interest, or that they also obey the master equations for the TMCs. Under such 
conditions the unfolding procedure then allows both the Fj^^\x) and the target mass 
contributions to be determined directly from data. The procedure is to parameterize 
fI;'^\x) at fixed Q^, insert this into the master equations, and then minimize the 
difference of the resulting F™^ with respect to the data. 

This method is complementary to PDF fits. One might expect the F2^^ structure 
function so extracted to be consistent with that obtained from PDF fits to the data 
when including TMCs. However, an advantage of the unfolding procedure is that the 
separation of the Fg^'' (which is described by pQCD) from the target mass contributions 
is good to all orders in a^, since Eq. (23) is valid to all orders. Comparisons of the results 
from the two techniques could help resolve any possible mixing between contributions 
from TMCs, higher twists, and those from higher order pQCD terms. Such an unfolding 
procedure has been undertaken [66] for the world data set of charged lepton scattering 
data from the proton, and the results are currently being prepared for publication. 

In this study the F2 data were fitted globally for 0.5 < Q'^ < 250 GeV^ by 
allowing for a dependence of the parameters describing F2^\x). The results for 
= 3 (top panel) and 20 GeV^ (bottom panel) are shown in Fig. 10. The solid curve 
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Figure 9. Comparison of the F2 structure function, with and without target mass 
corrections, and NuTeV data [64]. The base PDF set is CTEQ6HQ [7]. 



is the ^2°^ determine from the fit, while the dashed curve is the full F™"^. Consistent 
with the determination from PDF fits previously discussed, the TMC contributions to 
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Figure 10. Results of the unfolding procedure of Ref. [66] for the proton structure 
function, at = 3 GeV^ (top panel) and = 20 GeV^ (bottom panel) . The massless 
limit functions arc shown by the solid curves, and the target mass corrected results by 
the dashed curves. 



F2 are large at small Q^, as much as 9% even at X = 0.4 for = 3 QeVl While the 
TMCs become much smaller at higher Q^, they are still sizable at higher CC ^ clS CcLIl be 
seen in the inset in Fig. 10 (bottom panel). At = 20 GeV^ the contributions from 
TMCs are 4%, 8% and 14% at x = 0.65, 0.70 and 0.75, respectively. 

It is interesting to note that even in kinematic regions where the TMCs are large, 
the unfolding procedure gives results for the target mass contributions which are quite 
consistent with that determined by inserting existing CTEQ6 PDFs into the master 
equation for F2. This can be observed in Fig. 11, where the ratio of TMC factors 
determined from the unfolding procedure to those calculated from CTEQ6 PDFs are 
plotted. The correction factors are here defined as Fg^"* / F™'^, e.g. the multiplicative 
factor which would be applied to data to yield a measured F2'\ Figure 11 shows the 
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Figure 11. Ratio of target mass corrections for F2 determined from the unfolding 
procedure (global fit) of Ref. [66] to the corrections obtained from CTEQ6 PDFs cf. 
Eq. (64). The values shown are in GeV^. The x region for which > 4 GeV^ 
is indicated by the symbols. Note that even at high x and low the difference in the 
TMCs is typically less than several percent. 



double ratio 

) ' /"'-^"''^ (64) 

y z / z /cTEQe 

as a function of x for 3 < < 9 GeV^. For the region W"^ > 4 GeV^ (indicated 
by the symbols) the difference is typically less than a few percent even up to x ~ 0.9. 

As a final note, the size of the non-leading terms from the unfolding procedure is 
shown in Fig. 12 as a function of x for in the range 1 to 12 GeVl We display the 
ratio of the full F™'^ to that including only the leading term (top panel), and next- 
to- leading term (bottom panel); note the change of scale between the two panels. The 
fractional contribution from the leading term is found to be in excellent agreement to 
that shown in Fig. 8, while the contribution from the double integral term is at most 
1.5% at = 1 GeV^, and neglible for higher Q'^. 
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Figure 12. Ratio of the full F]-^'-^ , as extracted from the unfolding procedure, to 
that including only the leading term ^^mc, Leading ^^^^ panel), and next-to-leading 
term ^TMC,Ncxt-Loading (-i^q^^qq-^ panel), for = {1,5,12} GeV^. The x region for 
which W'^ > 4 GcV^ is indicated by the symbols. Note the change of scale between 
the two panels. 



6.5. Longitudinal Structure Function 

Finally, we show a set of plots for the longitudinal structure function, F^, in Fig. 13. 
The target mass corrected structure function F™'~'{x, Q"^) has been computed according 
to Eq. (26) using parton model results for the structure function F^\x., Q"^) for charged 
current neutrino-proton scattering in NLO QCD including quark mass effects. Because 
the massless leading-order parton model contributions to Ff^ vanish, this quantity is 
ideally suited to study both mass effects and higher order corrections. 

The effects of the TMC are illustrated in Figure 13a) where the ratio F^/F™^ 
is shown for = 1,4, 10, and 20 GeV^. Note that the present ratio is the inverse of 
similar ratios shown in Fig. 7 for the structure functions F2 and F3. As can be seen, 
the curves are steeply falling, implying very large target mass corrections. While the 
TMC are proportionately large in this region, unfortunately the longitudinal structure 
function becomes vanishingly small. This point is evident in Fig. 13b) where we plot the 
absolute magnitude of Fl. Here, the solid lines represent the full target mass corrected 
F™'^{x, Q"^), and the dashed lines represents the conventional parton mo 
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Figure 13. a) Ratio of the structure function Fl without (a;,Q^)) and with 
target mass corrections at = 1,4,10, and 20 GeV^ vs. x b) 
Absolute magnitude of the target mass corrected F]^^'~^{x, Q^ ) (solid) as compared with 
the structure function in the limit of a vanishing target mass Ff\x^Q'^) (dashed) at 
g2 = 4 and 20 GeV^ vs. x. 



We do observe for Q"^ = A GeV^ that there are significant differences in the intermediate 
X region (~ 0-3); however, these effects are neghgible for = 20 GeV^. In conclusion, 
we note that while the longitudinal structure function, Fi, is difficult to measure because 
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of its small size, a precise measurement of this quantity could prove to be sensitive to 
the TMC and quark mass effects. [67] 

7. Conclusions 

We have presented in this review a survey of the key issues pertaining to target mass 
corrections (TMCs) in inclusive lepton-nuclear scattering structure functions and their 
impact on the analysis of experimental results. As illustrated by our structure function 
results with and without TMC terms (Figs. 9-12), current experimental accuracy 
demands that these effects be properly incorporated. 

In the context of the operator product expansion (OPE), we have outlined the 
derivation of the master equation for the target mass corrected structure functions 
F™'^{x,Q'^), Eq. (21), without any assumptions about fixed-order perturbation theory 
or the Callan-Gross relation. However, when relating the OPE results to calculable 
parton model expressions, we are forced to make compromises; for example, we generally 
work with only at leading twist {of. Sec. 3) and at a finite order of as in perturbation 
theory. Because the TMCs enter as (M^/Q^)-' corrections, it is essential to organize our 
expansion so that these terms can be distinguished from true higher twist contributions, 
such as from four-quark operators. 

To illustrate the dependence on both hadron and parton masses in the formalism, 
we present explicit results for heavy quark production via neutrino-nucleon scattering. 
For the non-leading terms, we compute an upper limit for the integral terms, derive 
approximate expressions, and evaluate the full expression numerically; as expected, the 
dominant TMC effects arise at large x and small values. 

In the limit x ^ 1 we encounter an additional complication in the standard TMC 
formulation, associated with the non-vanishing of the target mass corrected structure 
function at x = 1. This problem arises if the massless limit structure function, 
-^i'^H^iQ^)) is not required to vanish in the kinematically forbidden region ^ > ^o- 
One may hope that higher twist terms could cancel the unphysical contributions. 
Alternatively, one can enforce the constraint that F2^\^,Q'^) does vanish for ^ > ,^0! 
which removes the unphysical contributions, but at the expense of introducing additional 
^0 dependence into the structure function. We have reviewed various attempts that have 
been made in the literature to address this problem. 

Finally, the importance of the TMCs for comparison with experimental data is 
established in Sec. 6. We perform fits of structure functions both with and without 
the TMC contributions, and observe that these terms can have a substantive impact 
in the region of high x. While there are a number of factors which enter these fits, 
it is encouraging to see that the TMC contributions generally improve the agreement 
between the fits and data. 

We also demonstrate (Figs. 10, 11, and 12) quantitative agreement extracting the 
F2 structure function with two complementary methods. We find that the ratio of the 
TMCs for F2 obtained starting with PDFs and including TMCs, vs. starting with the 
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TMC master equation, agrees at the 5% or better level, even for = 3 GeV^ for 
X < 0.65. 

The fact that the experimental data are sufficiently accurate to be sensitive to the 
TMC effects represents an important milestone in the study of inelastic lepton-nucleon 
scattering. From this foundation, future studies can more fully characterize the TMC 
effects, and begin to separately identify higher twist contributions for these processes. 
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Appendix A. 

Kinematics and Fundamental Relations 



Appendix A.l. 
Notation 

We consider the basic inclusive lepton-nuclear scattering process l{k) + N{P) 
i'{k') + X{Px) {cf. Fig. 1), with q = k — k' the four-momentum transferred to the 
nucleon, and define x to be the standard DIS scaling variable in the Bjorken limit: 

x = = , A.l 

2P-q 2Miy ' ^ ' 

where —Q^ = is the virtuality of the exchanged boson, and P is the target hadron 
four-momentum, with M the target nucleon mass 

(P2 = M^). We denote by nii the 
mass of the i'th quark. For a nuclear target with atomic number A, it is convenient 
to rescale the Bjorken variable x by A so that the denominator of Eq. (A.l) is still the 
mass of a nucleon. Additionally, we introduce the energy transferred to the hadronic 
system in the target rest frame as: 

iy = E-E' = ^ (A.2) 



and the inelasticity 

E-E' P-q 
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where E and E' represent the initial and final lepton energies, respectively. The mass 
of the hadronic final state is given by: 

iy2^P| = (P + g)2 = M2 + ^(l-x) . (A.4) 

For convenience, we also introduce the variable r, which is a combination of factors that 
appears frequently: 




Appendix A. 2. 

Generalized Nachtmann Variable 

For massless quarks, the parton light-cone fraction is given by the Nachtmann variable ^. 
Given the 4- vector = {Pt, Px, Py, Pz}, we can re-express this in light-cone coordinates 
as P^ = {P+,7^^,P_} wheretl P± = {Pt ± P.)/V2 and 7^^ = {P^^Py} are the 
(boost-invariant) transverse components. Light-cone coordinates are a convenient 
representation when the momentum components are strongly ordered, e.g. P+ ^ 
P±, P_. Thus, if the hadron light-cone vector is P^ = {P"*", , M^/ (2P"'")}, the parton 
vector for a (massless) parton is = {^P"*", ,0}.tt In the case of massive partons, 
the Nachtmann variable ^ is generalized to ^. [29, 57, 68] 

We identify two generalized Nachtmann-type of variables, ^ and ^, which are related 
to the Bjorken x variable via a dimensionless multiplicative factor. These can be written 
as: 

^ = xRm, (A.5) 

^ = ^R,,=XRMR^J , (A.6) 



where 



Rm = , = y— , (A.7) 

1 + Jl + (4x2MVg2) 1+r 



_ Q"^ - + + A(-g^ ml m^) 

R^J - ^2 ' 



A(a, b, c) = + 52 + ^2 - 2{ab + bc + ca) , (A.9) 

where mj is the initial quark mass and ruj is the final quark mass. The variable ^ = x Rm 
essentially corrects Bjorken x for the effects of the hadronic mass. Computationally, this 
arises from the final state momentum conserving delta-function, 5^(g + P — Px), which 
can be re-expressed to include the delta-function 6{x — ^/Rm)- 

In a similar manner, the variable ^ = ^ Rij further corrects the ^ variable for the 
effects of the partonic masses. The origin of this correction is the momentum conserving 

tt There are multiple conventions here; sometimes the 1/V2 is replaced by a 1/2 or a 1. 

ft The parton model is derived in the coUinear limit where the parton transverse momentum is neglected, 

e.g. ^ = 0; finite transverse momentum can contribute to the TMCs, c/. Sec. 1. 
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delta-function at the partonic level, S^{q + Pi — Pj), which can be re-expressed to include 

It is important to note that the hadronic correction Rm and the partonic mass 
corrections Rij are separate and can be applied multiplicatively. Also note that the Rij 
correction factor depends on the specific partonic masses involved; hence, this cannot 
be simply factored out as with the Rm terms. 

In the limit the initial quark mass vanishes (mj 0), the partonic correction factor 
has the hmit Rij ^ (1 + mj/Q^), which we recognize as the original slow-rescaling 
correction: 

Our focus in including quark mass corrections is for charm production in neutrino 
scattering, where incident quark masses are negligible, or where heavy quark PDFs 
are very small and additionally suppressed by small mixing angles. Consequently we 
shall focus on the mj = results below. 



Appendix B. 

Charm Mass Dependence in /i2 and g2 

We discuss here the charm mass, m^, dependence in the target mass corrected structure 
functions in Eqs. (44)-(46). Inserting Fg^"* from Eq. (40) into the definition of /i2(0 
Eq. (18), we find (for a diagonal CKM matrix): 

h2iO = / du ^ ; ^ = / du ^ = 2(1 + ml/Q^) du^ , 

u^ Je, M H u 

(B.l) 

where the variable u = u{l + ml/Q'^). Making a variable transformation to u, we find: 
h^iO = 2(1 + ml/Q') f'du'^ = 2(1 + ml/Q') du '^""^ 



U J£, U 

(B.2) 

Since the strange quark PDF vanishes for n > 1 we can restrict the upper integration 
bound to 1 and arrive at the final result: 

h^iO = 2(1 + ml/Q^) C du ^ . (B.3) 



u 



Similar steps lead to the results: 



92(0=2 tdu f'dv'-^, (B.4) 

h3{0= Cdu^. (B.5) 
Je u 
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Appendix C. 

Comparisons with the Literature 

In this appendix, we make the identification of our notation with the notation in earlier 
papers by Georgi & Politzer [34], De Riijula, Georgi & Politzer [40], and Barbieri et 
al. [35]. We also comment on the relation to the notation of the paper by Kretzer 
& Reno [42]. Known typographic errors are also indicated in this appendix, in an 
attempt to minimize confusion when comparing results in the literature. Since neutrino 
production of charm quarks has both target mass and charm quark corrections, we make 
the comparison for charm production in the diagonal Cabibbo matrix limit. 



Appendix C.l. 

Comparisons with Ceorgi & Politzer 

The results for the target mass corrected structure functions in the work of Georgi & 
Politzer [34] are given in their Eqs. (4.19), (4.20) and (4.22). Only leading order in QCD 
is considered, which in the massless limit yields the Callan-Gross relation for the Fi 
and F2 structure functions. As a result, a function -F(^) is the only function introduced. 
Before proceeding, we should correct two typographic errors in Ref. [34]: 

(i) In relation (4.19) of Ref. [34]: 

^'^^ '^^ " 2(l + 4MVg2)i/2^^^) + - ^ 2(l + 4a;2MVg2)i/2^(^) 

(C.l) 

(ii) In relation (4.22) of Ref. [34]: No derivation of W3 is provided to get the overall 
factor, but there is a factor 2 mismatch in relative normalization between the first 
and second terms: 

'^^3(Q',x) _ X X 

M 2(1 + 4x2 MVg2) (l + 4x2MVg2) 

(C.2) 

With these corrections, Eqs. (4.19) and (4.20) of Ref. [34] agree with Eqs. (22) and 
(23) at leading order if F{^) in [34] is identified with the structure functions F'f'\^) 
such that: 

F{0 = Fi'\i)/e = 2F?\0/i. (C.3) 

Note that F!f\^)/^'^ = 2Fi'^\c,)/C, is the Callan-Gross relation which holds in leading 
order for massless quarks. 

The target mass corrected structure function uW^/M found in Eq. (4.22) in 
Ref. [34], after correcting the typographic error, agrees with Eq. (24) if: 

= F(0 = 2Ff(0/e. (C.4) 

In Ref. [40] De Rujula et al. discuss the electromagnetic structure functions, where 
F{^) —>■ -Ff (0 or -^2^(0- The Ff are not defined, however, in terms of the structure 
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functions -F"/""*. In addition, the structure functions Wi and W2 are a factor of 1/M 
times the Wi defined in Ref. [34]. 



Appendix C.2. 

Comparisons with Barhieri, Ellis, Gaillard, & Ross 

In comparing with the notation of Barbieri, Ellis, Gaillard, & Ross [35], we note firstly 
that the quantity 'V in [35] is scaled by a factor M compared to our notation, and we 
label it with a subscript ub'- 

UB = iyM , (C.5) 

which gives x = Q'^/{2i/b) = Q'^ /{2Mu). Note also that the variable denoted by in 
Ref. [35] corresponds to the generalized Nachtmann variable ^ used in this paper. 
Substituting our notation in Eqs. (2.12a-e) of Ref. [35] yields: 

X 



B 



1 



B 



gv 



w^ 



Wi 



r^Q^' 



'IV 



2ab 



Ax^M^ 
r^Q^ 



'IV 



2xM^ 



B 



(C.6) 
(C.7) 



du 1 



Here = [{Q^ + m) 



■m 



2\2 



(C.8) 

-4m^(5^], rui and ruf are the initial and final quark masses, 
and a and h are the vector and axial-vector couplings, respectively. The structure 
function Ti corresponds to our -F/°\ and the Qi term contains the single and double 
integral terms corresponding to our and (72 expressions of Eq. (21). 



With mf <^ mj, in the first two structure functions 



(C.9) 



yields W^ and W2 with the same normalization as in Eq. (10). In the third structure 
function, we find that: 

s{i) 



2abJ^tiO 



(C.IO) 



We obtain agreement for the third structure function if ubW^ = 2M'^F™'^ . 

We note, however, that we do not match the expressions for Barbieri et al. [35] in 
the limit of a finite initial quark mass (m^ 9^ 0). While the initial quark mass terms are 
small numerically, for completeness we comment on the relation of our results to this 
work. 

The renormalization and factorization prescriptions allow certain freedoms in 
organizing the terms in the master equation; hence, the mass dependence need not be 
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unique. Our results are constructed to a) reduce to standard parton model expressions in 
various limiting cases {e.g. {mi,mf, M} 0), and to b) maintain the simple structure 
of Eq. (21). Our formulation incorporating quark masses follows the factorization 
construction of Collins [69]. In particular, the {A,B,C} coefficients for our master 
formula, Eq. (21), are free of any quark mass terms; the dependence is entirely 
contained in Pj^^K 

Our expressions also satisfy the expected limiting cases. For example, in the limit 
M — > our expressions reduce to the general parton model result with finite {mj,m/} 
terms, by construction; this is in contrast to Eqs. (C.6)-(C.8). 

In summary, our expressions of Eq. (21) follow the simple factorized structure of 
Ref. [69], and reduce to the standard parton model results in the limits of various 
vanishing masses. The numerical difference between this formulation and that of 
Ref. [35] is small. 



Appendix C.3. 

Comparisons with Kretzer & Reno 

In Kretzer and Reno [42], the notation has the slow rescaling variable including the 
charm quark mass written as ^/A, where A = (1 + m^/Q^)"^, so that ^/A = ^. The 
quantity p in [42] equals r, and /i = M"^ jQ^ . In the absence of charm mass corrections, 
the results of Ref. [42] are in full agreement with the equations in this paper. 

The translation of the treatment of _P^'^mc terms of Pf^^ in Ref. [42] is not made 
explicitly when charm production is involved. Charm production is discussed in terms 
of Ti ~ s. The expression for charm production in Eq. (3.27) of Ref. [42] is written as 
the sum of three terms times coefficients. The terms depend on jFj, TYj and Q^- The 
quantity J^j(0 appearing in Eq. (3.27) should read Ti{^, with 

m) = s{o (c.ii) 

at leading order. The functions Tij and Q2 in Ref. [42] for charm production should also 
be functions of ^ rather than ^. In terms of g2 and hi introduced here, one has: 

h2{0 = 2{l + ml/Q')n2{0 , (C.12) 
92(0 = 2^2(0 , (C.13) 
hsiO = m,{0 . (C.14) 
With these corrections, Eq. (3.27) in Ref. [42] is: 

The sum is from i = 1, ■ ■ ■ , 5, however, there is no mixing of the 4th and 5th terms in 
F™^. Tables IV and V with the coefficients a* and are correct, however, there is an 
error in Table VI. The coefficients 7^ all have an extra power of A in the denominator. 
Explicitly, the first two coefficients should read (with the translated notation): 

71 = ^ . (C.16) 
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7. = • (C.17) 

In fact, Eq. (C.15) was used in the numerical work in Ref. [42]. The contributions of 71 
and 72 are small, so the error in Table VI is not numerically significant. 
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